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Within the replica approach to mean-field spin-glasses the transition from ergodic high- 
temperature behaviour to the glassy low-temperature phase is marked by the instability of the 
replica-symmetric saddle-point. For general spin-glass models with non-Gaussian field distributions 
the corresponding Hessian is a 2" x 2" matrix with the number n of replicas tending to zero even- 
tually. We block-diagonalize this Hessian matrix using representation theory of the permutation 
group and identify the blocks related to the spin-glass susceptibility. Performing the limit n — > 
within these blocks we derive expressions for the de Almeida-Thouless line of general spin-glass 
models. Specifying these expressions to the cases of the Sherrington-Kirkpatrick, Viana-Bray, and 
the Levy spin glass respectively we obtain results in agreement with previous findings using the 
^ , cavity approach. 
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I. INTRODUCTION 



Spin glasses are paradigmatic examples for systems with competing interactions Both their equilibrium and 
dynamical behavior shows unique characteristics which are absent in systems without frustration. The concepts 
and techniques introduced in the theoretical description of spin glasses [2J have found interesting and widespread 
applications in other, at first sight unrelated fields of science such as complex optimization, error-correcting codes, 
artificial neural networks, and computational complexity 
. One of the central features of spin glasses is their non-ergodic low temperature phase characterized by slow relaxation 
and hysteretic response to external magnetic fields. A thorough theoretical understanding of this phase is available 
only for mean-field systems where the spin-glass phase is composed of a hierarchy of ergodic components. In the 
i parameter plane spanned by temperature and external magnetic field the high-temperature phase is separated from 

the glassy low-temperature phase by the so-called de Almeida-Thouless (AT) line [l], Q . The determination of the 
t-H ■ AT-line is therefore of central importance in the theoretical analysis of spin-glass models. 

Two rather different approaches are by now available to calculate the equilibrium properties of mean-field spin 
glasses. The replica method p] starts with n replicas of the system under consideration which after the ensemble 
average over the quenched disorder interact with each other. The free energy can be determined from a saddle-point 
integral over order parameters. The trademark of the replica method is the mathematically problematic limit n — > 
to be performed at the end. In this framework the AT-line is determined by the local stability of the replica-symmetric 
saddle point [f| ■ In the cavity method @, H| one spin is added to a system of N spins and the stochastic stability of the 
thermodynamic limit N — > oo is used to derive self-consistent equations for the order parameters. Here the AT-line 
may be obtained by investigating the correlations between two spins which must vanish in the thermodynamic limit 
• • . for a pure state of a mean-field system [2] . 

Both methods have been implemented for the analysis of the simplest mean-field spin glass, the Sherrington- 
Kirkpatrick (SK) model @. For this model the ergodic phase is characterized by a single order parameter and 
a Gaussian distribution of local magnetic fields. The fluctuations around the replica-symmetric saddle-point are 
described by an n(n— l)/2x n(n— l)/2 matrix. Its eigenvalues have been determined in [a, [13, 13 • The temperature 
dependence of these eigenvalues shows that the replica-symmetric saddle-point loses its stability at thephase boundary 
of the ergodic phase. The detailed form of the AT-line was reproduced within the cavity approach Q. 

The situation is less clear for more general mean-field spin-glass models which unlike the SK-model are characterized 
by non-Gaussian distributions of local fields. Models of this type are in particular important in complex optimization 
@,[I3- A prototype of this class is the Viana-Bray (VB) [l3[ model for a diluted spin glass in which each spin interacts 
with just a few, randomly selected other spins. Here the AT-line was determined numerically in (l4| . whereas analytical 
information is available only near the freezing temperature [l3j]. The replica treatment of diluted spin glasses and 
optimization problems is more complicated than that of the SK-model and involves already at the replica-symmetric 
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level an infinite number of ord er p arameters (la . 1 1 €31 ] . A general and elegant approach to this more complicated setting 
was introduced by Monasson [17j | . The fluctuations around the replica-symmetric saddle-point are now characterized 
by an 2™ x 2 n matrix which has to be diagonalized in order to assess the stability of replica symmetry. Recently 
it has been shown [18[ that this method may also be used to analyze spin-glass models characterized by coupling 
distributions with diverging moments such as Levy glasses [l9l l20j . This opens up the possibility to determine the 
AT-line also for such models within the replica method. 

In the present paper we investigate the stability of the replica-symmetric saddle-point for spin-glass models with 
non-Gaussian field distribution. To this end we implement the approach of Monasson for diluted spin glasses and 
reduce the determination of the free energy per spin to a saddle-point integral over 2™ order parameters. The Hessian 
matrix describing the fluctuations around this saddle-point can be block-diagonalized by exploiting the representation 
theory of the permutation group [2l[ . We also build on techniques introduced in [2^, for the analysis of replica 
symmetry breaking in one-dimensional spin glasses. We then identify the blocks which are related to the spin-glass 
susceptibility xsg the divergence of which signals the onset of spin-glass order. Up to this point the analysis is rather 
general and uses only the replica structure of the fluctuation matrix. The final diagonalization of the relevant blocks 
can only be performed after the details of the model under consideration are fixed. We consider three representative 
examples: the SK model which merely serves as test case for our method, the VB model as example for diluted spin 
glasses, and the Levy glass as system with a local field distribution exhibiting long tails. In all cases we provide 
expressions for the AT-line separating the replica-symmetric part of the phase space from the region characterized by 
replica symmetry breaking. 

The paper is organized as follows. In section [TT| we define the central models of interest, recollect the main steps 
in the replica-symmetric theory for diluted spin glasses and fix the notation. Section IIIII contains the analysis of the 
situation without external magnetic field for which the calculations are significantly simpler. Section IIVI is devoted 
to the general case from which the expressions for the complete AT- lines in the models considered result. Finally, 
in section [V] we give some conclusion and discuss open problems. Some more technical steps are relegated to the 
appendices. 



II. BASIC EQUATIONS 

We consider Ising spins Si = ±l,i = 1, N with random, pairwise interactions specified by a symmetric matrix 
Jij in an external field /i ex t ■ The Hamiltonian is of the general form 

H ({Si}) = ~\Ys J :r S ' S > ~ h ^ E S * > W 

(*>i) ' 

where the first sum runs over all pairs of spins. The couplings Jij = Jji are i.i.d. random variables drawn from a 
model dependent distribution Pm- We will consider three examples for this distribution in detail, namely 

fW ( NJ 2 \ 

Psk(J) ■= Jf- exp — ) (2) 



47T \ 4 

Pvb(J) ■= ^p(J)-(l-^)s(J) (3) 

£,(?):= J dJe- u "P a (J) = exp(-^j a 6(0, 2] (4) 

corresponding to the Sherrington-Kirkpatrick (SK) model [9], the Viana-Bray (VB) model [13j . and the Levy spin 
glass [18[ respectively. Here p(J) denotes the distribution of the non-zero bonds in the VB model. The distribution 
P a for the Levy spin glass is defined via its characteristic function P a . Note also that the variance of the Gaussian 
Psk is twice its standard value such that it coincides with the Levy case in the limit a — > 2. 

The SK model is the paradigmatic case of a fully connected spin-glass in which each spin interacts with each other 
via weak couplings of order N~ x / 2 . Complementary, the VB model is characteristic for diluted spin glasses in which 
each spin interacts with only a few other spins via strong 0(1) couplings. The Levy spin glass interpolates between 
these two extremes since each spin interacts with each other spin but most of the couplings are very weak whereas 
0(1) couplings per spin are strong. 

The large connectivity limit of the VB model leads either to SK- or Levy-like behaviour in the thermodynamic 
limit, depending on the existence of the second moment of the distribution p(J). If the second moment exists this 
limit is defined as 

k i y iV, p{J) i y y/Np(JNJ\ , (5) 



3 



and leads to an SK-modcl due to the central limit theorem. On the other hand, if the distribution p(J) has a power 
law tail: p(J) — |j|^oo I J|~ ( - Q+1 \ a < 2, the large connectivity limit is given by: 

k^N, p(J)^N°p(n°J\. (6) 

The limits coincide for a = 2 as it should be. 

In the framework of the replica approach the free energy density /, is expressed through the disorder averaged 
replicated partition function Z n (j3) via 



1 Zn (P) 1 



f:=- lim — InZ (/3) = - lim — lim ^ , (7) 

where for integer values of n the replicated partition function is given by 

Z n 03) = ]T exp (-/3 H ({St})) . (8) 

After performing the disorder average in ((5|) the trace over the spin variables can be transformed into a 2™-dimensional 
integral over order parameters [l7T | 

c (*) = "i!?E ( 9 ) 



N ■ 

where a = {<J a } a =i.„ n stands for an Ising spin vector with n components and 5(S,<r) denotes the Kronecker-£. Hence 

(10) 



(Z n {P)) j = J n dc (3) 6 c (3) - l^j exp (-Npf trial ({c (a)})) , 



where the 5-function enforces the constraint 

X>(<r) = l (11) 

resulting from (jH]). The trial free energy f trial has two contributions f trial = fs + /s according to 
/3/tria! ({c(<?)}) = c (ct) In c (ct) - 



i^c(a)c(f) J dJG M (J) exp((3Ja-f)+f3h cxt ^c(ct) (X^) 



The first term is the entropic contribution fg specifying the number of spin configurations realizing a particular 
set of order parameters c(B). The second term /e derives from the Hamiltonian ([1]) and comprises the interaction 
energy and the energy in the external magnetic field. This latter contribution as well as fs do not depend on the 
explicit model considered. The interaction part in f& on the other hand depends on the specific form of the coupling 
distribution Pm which is encoded in Gm(J)- For the three models specified in ©-((I]) one has 

Gsk(J) =S"(J) Gvb{J) = k(p{J)-5{J)) G l {J)=-J ^ | J\ a exp( i J J), (12) 

respectively, where 5" denotes the second derivative of the (^-function. In the case of the Levy spin glass we work at 
imaginary temperatures j3 = —ik as long as n =/= 18]. 

With the form (flO|) the calculation of the free energy / is reduced to a single site problem as is characteristic for 
mean-field systems. In the thermodynamic limit N — ¥ oo the remaining integrals over the order parameters can be 
evaluated by the saddle-point method. One therefore has to determine the minima Co of f trial ({c(B)}) satisfying 



cp (g) = £ n exp f cp (r) J dJ G M (J) exp(/3 JB-f) + ft ftext Q"q J , (13) 

where C n is a factor accounting for the constraint (|llj) . 

The general solution of the saddle-point equation is a complicated problem. To get some guidance one uses as a 
first step the assumption of replica symmetry (RS) stipulating that cq depends on B only through the sum J3a=i aa °f 
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the vector components. These replica-symmetric order parameters are related to the (replica-symmetric) distribution 
of local fields V(h) in the spin glass by [17] 

As is well known the replica-symmetric solution fails at low temperatures which on the formal level is due to the 
instability of the replica-symmetric saddle-point. To assess this stability the temperature dependent eigenvalues of 
the Hessian % describing the quadratic fluctuations around the replica-symmetric saddle-point have to be determined. 
The matrix elements of T-L are given by 



r-ti / —> — #\ P d 2 ftrial 

n{<T > T) = dc(a)dc(r) 



% if --r- fdJGu(J) exp(^Ja-r). (15) 

0=1 a a) J 



.EL 



Similarly to f trial the Hessian splits into an entropic and an energetic contribution, H = Hs + He- Note that the 
external field does not show up explicitly in H. Nevertheless it influences the stability of the replica-symmetric 
saddle-point due to the dependence of Co on h cyLt as specified by (jT5)) . Note also that the expression (|T5|) does not yet 
account for the constraint (fTT|) . 

The determination of the eigenvalues of % may be systematically simplified by exploiting the permutation symmetry 
of the replica-symmetric saddle-point. To this end it is convenient to express the Hessian as a tensor product of 2 x 2 
matrices using 

= f[ =f[(a a \ (e^«7 + e"^) \r a ) = (a\ (g) (e^a + e^ J a x ) |f) (16) 

a— 1 a— 1 a— 1 

for the energetic and 

-iJ^-^f ^ e ^f[^e-^={t\ I ^^(g)expH^ 3 )|r) (17) 
c (Ea=i<M J 2tt c(r) £± J 2tt c(r) ^ 

for the entropic part, respectively. Here the vectors \a) are defined as \a) — \a±, oi, . . . <r„) — \ <J a ) with a a = ±1 

which span the space V of replicated spin configurations. The matrices 

(18) 





















<?0 = 


i 






! 


1 


, ct 3 = 


i 





are the usual Pauli matrices. 



III. ZERO EXTERNAL FIELD 



The discussion of the eigenvalues of the Hessian T-L is relatively straightforward for /i ext = and T > T c . We 
therefore first derive expressions for the eigenvalues in this region and afterwards turn to the more involved situation 
with /i oxt ^ 0. 

For h ext = the saddle-point equation (fT5)) has for all j3 the solution 



(19) 



which, according to HU, corresponds to P(h) — 5(h). On physical grounds we expect that this paramagnetic solution 
is stable at sufficiently high temperatures. For constant c(s) the Hessian (fT5|) for the paramagnetic solution reads 

n r. n 

(a\ Hpm \f) = (a\ 2" (g) a - / d J G(J) (g) (e^ J a + e^ai) |r) . (20) 

ct=l ^ a=l 

Its eigenvectors can therefore be constructed from those of cti, namely 

|fc>-^(l+} + (-i) fc |->) fee {0,1}. (21) 
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The eigenvectors of Upm may hence be written in the form 

n 

\k) = \h,...,k n ) =®|*„> with |A a ) G{|1),|0)} . (22) 

a=l 

The 2™ mutually orthogonal vectors \k) form a basis in the space V . With the abbreviation k — we find 

n 

(g) (e^&o + e'^ai) \k) = (2cosh(^J))"-' £ (2sinh(/3J))' £ (23) 

a=l 

and therefore the eigenvalue corresponding to |fc) is given by 

k l n = 2 n (l- [ dJG M (J) cosh"(/3J)tanh fc (/3J) ) . (24) 



Not all of these eigenvalues are, however, relevant for the stability of the saddle-point since the constraint (fTTj) needs 
still to be taken into account. To do so consider small fluctuations Sc(a) around the RS saddle-point cq{<j) = 2~ n . 
With the notations 8c(a) — {8c\a) and Sc(k) — (5c\k) we find from (|T0l) to second order in Sc(a) 



Z»(fi) « exp(-N(3f trial ({2- n })) Jl[ddc(a) si X>(<?) j exp I —^ScffiHpM&tfSct?) 

a \ & / \ o , t j 

= exp(-iV/3/ t rt /({2-"})) J l[dSc(k)S (2*<Jc(0)) exp j-y £ A^c(fc) 2 



(25) 



exp 



(-W trio/ ({2-"})-|ln2) | pd5c(fe) exp I £ A £< 5c(fc) 2 j 

fc#0 \ fc#0 / 



Here we have used the transformations 

*c (a) = £ (Sc\ a)=J2 <<M k) \S) = Sc(k) 2? % g = 2* £6(0) 

ct e? k ® k 

and 

£ 5c (a) r)5c (f) = £ Sc(kf . 

Consequently |0) is perpendicular to the constraint <jTTJ) and the integration in the corresponding direction is suppressed 
by the (^-function in (fTU| . The value of A° is therefore not relevant for the convergence of the integral (l25l) . 



(k) T. 

The limit n — > can be performed now for every eigenvalue A„ := A„ corresponding to the eigenspace spanned 
by eigenvectors with magnetization k 

K% = lim A* = 1 - / d J G M {J) tanh fc (f3J) . (26) 

The SG transition is signaled by the divergence of the SG susceptibility given by [l[ 

XSG := i £ <^>I lim — £ (StSfSlSS)^ . (27) 

*j *>i (o,6) 

where the second equality holds only in the paramagnetic phase. The last average is to be taken with respect to all 
distinct replicas after the disorder average has been performed and yields 

(S^4> repl = ]T SfS^exp I / dJG M (J)ef> jS *- 3 ' + Ph ext J2J2 S ?) . (28) 

{S«} \ *J i=lo=l / 
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Using the same method as for the replicated partition function Z n (j3) it is possible to rewrite ([28} as an 2 n -dimensional 
integral. In the saddle-point approximation one finds using cq(o) — 2~™ 



N 

XSG = N + hm — 

n->0 nyn — 



— J IJ dc (?) s( £ c (?) -l)j2c(a)c (r) (a ■ ff exp ( - N(if tnal {{c (a)}) 



— <^Tj — 7 II dSc W 6 [ E fc (*) J E Jc (*) • ex P hy E <K?)?W?, 



a .t 



1 

~ W)' 

The spin-glass susceptibility is hence directly related to one of the eigenvalues of the Hessian. As expected the 
divergence of xsg at the spin-glass transition corresponds to the instability of the replica-symmetric saddle-point 
signaled by A^ 2 -* = 0. We also mention that for non-symmetric coupling distributions a transition to a ferromagnetic 
phase may occur. This transition is related to the ferromagnetic susceptibility 

X F m := i £ WW. = ^ E £ E = X^) (29) 

and its divergence is hence linked with another eigenvalue of %pm- 



Examples 

We now discuss the results for the three particular models defined before. For the SK-model we find from ([26]) 
A^-l-2<S 2 , fe /3 2 . (30) 

In accordance with the symmetry Psk(-J) = Psk(J) (cf. ©) only a transition to a spin-glass phase is possible at 
low temperature. Defining the transition temperature T c by AS = we find T c — \f2 and 

T 2 



XSG 



T 2_ T 2 



which coincides with the well-known results 

The expression for the eigenvalues of the Hessian for the VB-model 

= 1 - k J dJp{J) tanh fc ((3 J) + n8 k . (31) 

was already derived by Monasson [2^]. Since tanh(/3J) < 1 the integrals become smaller with increasing k while 
the fluctuations corresponding to k = are suppressed by the constraint (jTTJ) . For symmetric distributions p(J) the 

(2) (2) 

integral for k = 1 vanishes identically and hence A^ is the first eigenvalue to become negative. From A) J, = results 



the well-known expression for the critical temperature of the spin-glass transition in the VB-model |13l [l5l fl(| : 
1 = k J dJp(J) tanh 2 (^ c J) (32) 
For the Levy glass we have 

Ait„ = l + / ^-\J\ a exp( l JJ)tanh fc (/3J) . (33) 
Again P a {J) = P a {—J) ensures that 

. (?) f dJdJ , ~,„ , ~. 9 . T(a + 1) /oar\ f dJ 9 , , 

Ail, = 1 + _/ I J| Q exp( i J J) tanh 2 (/?J) = 1 - ^ - ' sin (— ) J -— r tanh 2 (j5J) (34) 

is the first eigenvalue to become negative. For the critical temperature of this spin-glass transition we hence find 

T o = C(a) J tanh 2 (J) with C(a) = ^±12 sin ^ (35) 
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which was previously derived by both the cavity [19| and the replica method [18|. For the SG-susceptibility of this 
model results 

XSG — rp a _ rp a 

which is similar to the expression in the SK-model and in fact comprises it for a = 2. 

We therefore reproduce for all three examples the known results about the instability of the paramagnetic phase. 
However, it remains unclear at this point whether these instabilities can be cured within the replica-symmetric sector 
or whether replica symmetry breaking is necessary to stabilize the saddle-point. This question can only be addressed 
by investigating the more complicated situation with a non-trivial form of c( a a ) to which we turn now. 

IV. THE GENERAL CASE 

At low temperatures or /i ext ^ the Hessian % explicitly depends on a$ and the determination of its eigenvalues 
becomes much more complicated. The essential steps of the analysis are as follows: We first use the permutation 
symmetry between replica indices characteristic for any RS solution to block-diagonalize the Hessian according to the 
irreducible representations of the permutation group. For the remaining diagonalization inside the blocks the limit 
n — > can be performed and the eigenvalue equations assume the form of integral equations. We then show that 
the replicon eigenvalue related to the spin-glass susceptibility can only belong to one of the first three blocks. By an 
expansion around T c we then verify that all eigenvalues of the zeroth and first block which go to zero at T c return to 
positive values below T c . Consequently the replicon eigenvalue must lie in the second block. Finally we identify this 
eigenvalue and give explicit expressions for the AT-line for the three special models considered. 

A major simplification of the general eigenvalue problem of the Hessian is obtained by invoking the symmetry of H 
under permutations of the replica indices. Formally this symmetry is expressed by the commutation of the Hessian 
with a representation D of the permutation group S n . In our case D is defined by 

D(ir) \<T 1 ,<J 2 ■ ■ .(T n ) = 1^(1), 0^(2) • • -O-jr(n)) Wa) € { |+) , |-> } 7T G E n , 

which clearly commutes with % , i.e. 

D(tt)H = HD(tt) V7reS„. (36) 

The Hessian can therefore be block-diagonalized according to the irreducible representations D^ p > of E n that are 
contained in D. As a first step of the reduction of D we note that the subspaces V^ a > , spanned by vectors \ai, cr 2 . . . a n ) 
with a fixed number a of entries a a equal to plus one, are invariant under D. The restricted representation A^) of 
D on the subspace a = . . . n can be further decomposed into irreducible representations: 

AW~#e#-ffi])( s ' (37) 

with a = min{<7, n — a}. The irreducibihty of the D^'s is shown in [2l|. Each representation gives rise to an 
eigenvalue A' p ) of "H, with degeneracy 

deg (A<">) = dim (DW) = M - ( n _ \ , (38) 

The subspaces hence split into direct sums of irreducible subspaces V^'^ each of which is associated with a 
representation 

V W ~ 0) V {a ' 1} • • • © V {a ' d) . (39) 

This decomposition can be accomplished by Young-symmetrizers constructed using the Young tableaus 24]. The 
procedure shows that the vectors 

\\*;p)):=(\+)\-)) P \+r P \-) n - a - p e^ 1 , p = 0...a, (40) 

are mapped onto the irreducible invariant subspace V^'^ by an anti-symmetrization in two successive entries in the 
first 2p arguments, and a complete symmetrization in the last n — 2p entries. The vector 

\a;p) := A 2p S n - 2p (a) \\a; p)) = ( |+) R - R |+) ) P \a p) n _ 2p (41) 
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hence lies in V^ cr,p ', where the operators Av p and S n -2p(&) denote the anti-symmetrization and the symmetrization 
operators, respectively and the symmetrized part of the vector \a; p) with a — p entries equal to plus one is denoted 
by 

n—2p 

W~P)n-2 P = E ®K>- ( 42 ) 

^2<j a =2a-n a=l 

A basis of the subspace V^'^ can be constructed by applying all the D(tt) on \a; p) and choosing a maximal linearly 
independent subset. We note that the vectors {|cr; p) , a = . . . n, p = . . . a} are orthogonal, but not normalized. 
For a fixed p the set of the normalized vectors 

f |CT;P> cr = p...n-4 (43) 



V top 



is an orthonormal basis of a subspace of V, containing one element from each irreducible subspace 

V^ a ' p \ u = p . . . n — p. The matrix 'HgJ m with matrix elements 



< m (^r) = „ ====== : < r,r = p...n-p (44) 



/(p) c„ ^ _ (g; p| ~H |r; p) 

Vtopl o-;p) topi t;p) 

therefore comprises information from each irreducible subspace V < - a ' p ' > , a = p . . .n — p [22|, [2^, [25[ . Diagonalization 
of the n + 1 — 2p dimensional matrix Wg™, leads to the eigenvalues 

At p) , v = p...n-p (45) 

each of which is associated with one representation Z)( p ) arising in the decomposition of D. 

Exploiting the symmetry of the problem we have hence reduced the 2™ dimensional eigenvalue problem to I 5j 
eigenvalue equations of dimension n + 1 — 2p which are parametrized by p. Here [x\ denotes the largest integer smaller 
than x. As shown in appendix |B] the permutation symmetry of T-L can be further used to switch from the eigenvalue 
problem with symmetric matrices 'H^ Tn to non-symmetric matrices with matrix elements (a, t) = ((rr; p\\ H\r; p). 
This form allows an elegant continuation n — > which was already used in (2^, . We expect that the eigenvalues of 
the matrices T-L^ are generically non-degenerated since the permutation symmetry is already completely accounted 
for. 

Finally the eigenvalues corresponding to the representations have to be determined from the eigenvalue 
equation: 

n— p 

A (p) x (p) (r) = «W (a, r)x W (a) (46) 

a=p 

The decomposition T~L = Hs + He of the Hessian still holds in the new basis. The entropic part Hs is diagonal and 
depends on the RS solution c( a a ) of the saddle-point equation only: 

n { s P \a,r) = -^- (47) 
c(2er — n) 

The energetic part %^ depends on the details of the Hamiltonian and will be specified in the treatment of the 
particular models below. For the following analysis it is convenient to multiply the eigenvalue equations with the 
inverse of W£ and to transform the eigenvectors to functions <^ p ' according to [22|, l23j 

n — p 

<f>^ (h) := J2 X (P) W) exp (Ph(2a - n)) . (48) 

<T=p 

The limit n — > may then be performed which transforms the finite dimensional matrix equations into integral 
equations. In the case of the Levy glass it also allows the continuation to real temperature. Some intermediate steps 
of the calculations are given in appendix [CJ 

Eventually we arrive at eigenvalue equations of the form 

A (p) / AhV{h)(t) {p) {h + h) = cj) (p) {h) - [ dhV{h) [ d,JG M (J) {d h uf (h + h, J) cj) (p) (u(h + h, J)) (49) 
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where u(h, J) = -gatanh( tanh(/3/i) tanh(/3J)) and V(h) denotes as before the replica-symmetric distribution of local 
fields. This equation is the central result of the present section. 

As a first test we reproduce the spectrum of the Hessian obtained in section Hill for h cxt — and T > T c . In this 
case ^(h) = 8(h) and the eigenvalue equations simplify to 



\W<j>W(h) = <j) {p) {h)- j &JG M {J) {d h u) p (h,J)<p^{u{h,J)). (50) 
Setting h = and using u(h = 0, J) = as well as dhu(h = 0, J) = tanh (j3J) we get 

A<">0&O(O) = 0W(O)- / dJG M (J) tanh p (/3J)0^(O). (51) 



If ^W(O) 7^ we hence find A^ p ) = A^ p ^ and therefore reproduce expression (|26p for the eigenvalues obtained more 
directly in section IIIII If on the other hand 4>^ (0) = 0, then eq. (|5Tj) does not convey any information about A^ p -* . 
However, in this case we find after differentiating (1501) with respect to h and setting h = afterwards 

A (p) (0 (p) )'(O) = (^ (p) )'(0)- fdJGu(J) tanh (p+1) (/3J)(<// p >)'(0). (52) 



If (<p^) (0) 7^ this returns the expression for A' p+1 ). If (</> ( - p - ) ) (0) = we turn to the second derivative of (fBTj) and 
so on. In conclusion we find that for /i cxt = and T >T C the sub-block defined in (J3U) generates all eigenvalues 
y\W f rom (pgj) with k > p. Conversely the eigenvalue A' 2 ) which signals the spin-glass transition shows up only in 
blocks p — 0, 1, 2 implying that in some neighbourhood of T c all eigenvalues X^ with p > 2 are strictly positive. This 
is also corroborated by a replica representation of the spin-glass susceptibility starting with (|27[) which shows that 
the SG susceptibility does not depend on eigenvalues A^ p ^ with p > 2. In the following we therefore investigate only 
the "dangerous" blocks p = 0, 1, 2. 

For p = the constant function (fy-°'{h) = 1 is an eigenfunction corresponding to the eigenvalue A^ -* = 1 = A^ ). 
In the high temperature region this eigenvalue was irrelevant for the stability due to the constraint We assume 
that the same holds true in the spin-glass phase as well. For the first derivative of (|4"9")l we find for (p = 0) 

A (0) fdhV(h) (<f>W)'(h + h) = (<f> (0) )' (h) - [dhV(h) [ dJG M (J)d h u(h + h,J)((t>W)'(u(h + h,J)). (53) 



Hence either (<^°)) ee or is an eigenfunction of (|49)l for p = 1. Conversely if cf) 1 - 1 ' is an eigenfunction of (|49|) 

for p = 1 its primitive "I^ 1 ) satisfying (<I>( 1 ') / = fulfills the equation 

h {A (1 > / d^(/i)$ (1) (^ + /i)) = 9 h ($ (1) (^)- / d/*P(/i) / dJG A/ (J)$ (1) (u(/i + /i, J)) j . (54) 



Integration of this equation yields an eigenfunction <I>( 1 ) of (|4"9")) for p = since the integration constant may be 
absorbed in the choice of Consequently the block with p = contains the same eigenvalues as the block 

corresponding to p = 1 and in addition one eigenvalue corresponding to a constant eigenfunction which we believe to 
be irrelevant due to the constraint (fTTj) . This degeneracy between the p = and the p = 1 block is similar to the 
well-known degeneracy of the longitudinal and first transversal eigenvalue in the stability analysis of the SK model 

1- 

We now show that the eigenvalues and A^ which are degenerate with A^ 2 ) in the high temperature phase 
return to positive values below T c . In view of the equivalence between the eigenvalues from the zeroth and first block 
it is sufficient to show this for X^K We study the eigenvalue equation perturbatively to leading order in the reduced 
temperature r = 1 — T/T c at zero external field. To this end we expand the derivative of the p = 1 eigenvalue equation 
with respect to h at h — up to order h 2 [27| . The integral equation acquires the form of a 2-dimensional matrix 
eigenvalue problem. To leading order in t we find 

A«=2r, X^=t' 2 r, X^=ar, (55) 
efined as 

d /* , . u2 / p c J 



where t' 2 is defined as 



<2 := ~~ k J dJp(J) tanh^ [^f—^j > 0. (56) 

In all three cases the eigenvalue hence returns to positive values. The instability of the paramagnetic saddle-point 
due to unstable directions from the zeroth and first block are therefore cured by the replica symmetric low-temperature 
solution. The "dangerous" direction related to the replicon eigenvalue is contained in the p = 2 sector. Its detailed 
discussion requires a specification of the Hamiltonian which we therefore perform separately for the three case of 
interest. 
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The SK Model 

Using Gsk(J) = 5"{J) in 09]) for p = 2 wc find: 

J dhV(h)^(h + h) = ^(h)-2(3 2 ^(0) J dhV(h) — 4 ^ {h + ~ h)) ■ ( 5? ) 

At zero external field close to the transition temperature this equation can also be studied perturbatively in the 
reduced temperature r. An expansion of the last equation at h = up to order ft 4 turns the integral equation to a 
three dimensional eigenvalue problem. One eigenvalue becomes negative: 

A< 2 > =-^t 2 + 0(t 3 ), (58) 

indicating the well-known instability of the RS solution. 

In the presence of an external field we identify \^ = with the instability-line, which starts at T c for h cxt = 0. 
The replica symmetric distribution of local fields is explicitly known for the SK model: 

W^expf-fc^V (59) 



4irq \ Aq 

where in the spin-glass phase q is the non-zero solution of q = j dhV(h)tanh 2 ((3h). Setting A^ 2 ) = and h = in 
([57)) we arrive at 

1 = 2 B* eXp t~^ (60) 

J V2^ cosh 4 (P(V2qx + h cxt )) 

which reproduces the famous AT-line for the SK- model in the h ext -T-pla.ne 

The VB Model 

In the case of the VB model the eigenvalue equation for p — 2 reads: 

A (2) J &hV{h)(j> {2) (h + h) = <j)^(h)-K J dhV(h) J d,J p{J) {d h uf (h + h, J) <j>W (u(h + h, J)) . (61) 
At zero external field an expansion of the eigenvalue equation (1611) can be invoked leading to the eigenvalue 

A« = -I'-^^^ + Oir 3 ) (62) 

with 

t A := k / dJp{J) tanh 4 (^ c J) < k / dJ p(J) tanh 2 (^ c J) = 1 , 



and t' 2 defined in eq. (I5B1 . Due to the inequality in the last line is negative below T c indicating the instability of 
the replica symmetric solution for this model. 

Within the cavity approach the AT-line is described by [28[ 



fifth) = k J dhV(h) J dJp(J) (d h uY (h + h, J) cj)(u{h + h, J)) , (63) 

where the RS phase becomes unstable when the largest eigenvalue fi exceeds the value 1 [291 ] . Since (EU) for A< 2 ) = 
and (|63p for p — 1 coincide we have reproduced the stability criterion from the cavity method within the replica 
approach also for non-zero external field. 

The largest eigenvalue of eq. (|63[) can be determined numerically by simple iteration. Unfortunately we do not 
know about a similar straightforward method to determine the smallest eigenvalue A^ 2 ^ of eq. (|6Tj) . 
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The Levy spin glass 

The eigenvalue equation for the Levy SG is similar to the one for the VB-model: 

\W J &hV(h)(t> (2) (h + h) = <t> m {h)-C{a) J dhV(h) J (d h uf (h + h, J) (u(h + h, J)) , (64) 

where C(a) was defined in (|35|) . The expansion of the eigenvalue equation (j64")) to the leading order in the reduced 
temperature amounts to 

A(2) = ^1 + 2^ ^ (65) 

with 

U, a :=C(a) J ^^tanh^J) < C{a) J tanh 2 (/3 c J) = 1. 

This proves that A*- 2 -* is indeed negative below T c and the replica symmetric phase is unstable below T c . 
The stability analysis for the Levy SG performed in |28| using the cavity method gave rise to the equation 



fifth) = C(a) / dfiV(h) / rjr^Tt (dhuf (h + h, J) ftu(h + h, J)) 



\J\ 

where the instability of the RS solution was again signaled by fi > 1 30]. In the presence of an external field we 
therefore find equivalence between the results obtained using the cavity and the replica method. 

In the large connectivity limit of the VB model and in the SK-limit (a — »• 2) of the Levy SG all our results are 
consistent with each other. In view of (fl~2"|) we have lim Q _j.2 Gx(J) = 8" (J) — Gsk{J)- To obtain the large connectivity 
limit of the AT-line for the VB model we use © and ^ respectively in In the limit N —> oo the eigenvalue 

equation then acquires the desired form up to a constant depending on the details of the distribution p(J) which can 
be absorbed in the energy scale. 

To see the equivalence for the eigenvalues close to the transition temperature we use 

t' 2 ^a^2 and t A -> t 4 , Q -> t 4 , 2 = / dJ5"(J) tanh 4 (/3 c J) = (66) 

where the first arrow corresponds to the Levy limit of the VB model, and the second one to the SK-limit of the Levy 
SG. If the second moment of the distribution p(J) exists, one obtains using ([5]) directly t' 2 — > 2 and tn — > 0. 



V. SUMMARY 



In the present paper we derived within the replica formalism expressions for the AT-line of general mean-field 
spin-glasses including strongly diluted and Levy spin glasses. Due to the non-Gaussian character of the local field 
distribution in these models an infinite number of order parameters is needed already at the replica symmetric level. 
Following the approach of Monasson the fluctuations around the replica symmetric saddle-point are described by an 
2™ x 2" Hessian matrix. 

At high temperatures and in zero external field the distribution of local fields is a delta-function and the deter- 
mination of the eigenvalues of this Hessian is relatively straightforward. We find that all eigenvalues are positive at 
sufficiently high temperature and that some of them tend to zero at the critical temperature, T c , which signals the 
transition to the glass phase. 

Below the critical temperature the RS order parameter develops a non-trivial structure and the determination 
of the spectrum of the Hessian becomes rather involved. However, using the symmetry of the saddle-point under 
permutations of the replicas the Hessian can be block-diagonalized and the sub-blocks relevant for the stability of RS 
can be identified. Performing the n — > limit in these blocks turns the finite dimensional eigenvalue equations into 
integral equations from which general expressions for the AT-lines may be derived. 

We show the validity of our approach by applying it to three representative model systems: the SK-model as 
the standard model with Gaussian field distribution, the VB-model as the standard model for diluted spin-glasses for 
which higher moments of the field distribution are essential, and the Levy spin glass as standard model for spin-glasses 
with diverging second moment of the coupling distribution. 

We believe that with our stability analysis the replica-symmetric theory of general spin-glass models is now complete. 
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Appendix A: The spectrum of 7^' p ' 

First of all we show that the eigenvalues of n$y m and coincide. To this end we note that the symmetrization 
operators S n -2 P {cr) and Ai p can be represented as appropriate combinations of D(tt), tt G £„. The operator S n -2 P (<r), 
as a sum of ( n ~_^p) permutations, and A2 P as a product 

p 

A ?p = II [D(e)-D((2a,2a-1))] , (Al) 

a=l 

where e denotes the identity of the group and (2a, 2a — 1) the transposition of the elements 2a and 2a — 1. Being 
elements of the group algebra all operators S and A commute with H. Therefore 

{a; p\ H |r; p) = {a; p\\ H A 2p S n - 2p {a) |r; p) . (A2) 

The vector |t; p) is symmetric in the last n — 2p entries. The action of any D(tt) acting only on these last entries is 
hence trivial and we find 

S n . 2p (a)\r;p)= (™ ~ 2 j) h p) , (A3) 

since 5„_ 2p (cr) is a sum of (7_ 2 p ") permutations. From [D(e) - D((l, 2))] ( |+> |-) - 1-) |+) ) = 2( |+) |-> - |-> |+> ) , 
one similarly derives A2 P |r; p) = 2 P \r;p). The action of A2 P S n ~2 P {o) reproduces |r; p) up to the constant 2 p (™~ 2 p p ) 
which is the squared norm of \a; p). One hence has the following relation between the matrix elements of "H and T-L^ p ': 
{a; p\ H |t; p) = {a; p |cr; p) H^(<j, r). Using this relation in the eigenvalue equation for "Hg^ m : 

n-p n~ P r : ~ 

A W X W (T) = ^< m (a,r)^H = E\/^M^ (P) (^r)x (rf (-) 

O — p (T — p V \ 5 Pi 

& \^x ip) (r) = J2h^(<j,t) x ^(t), with x W(<7) = V(<7;/>k;p)x (p) W 

o= P 

we see that the eigenvalues of the matrices coincide. 



Appendix B: The limit n — > 

To perform the limit n — > we switch to the characteristic functions defined in (|48[) . Transforming the whole 
eigenvalue equation (l46l) one has to calculate the quantity 

n — P 

^H (p) (CT,r)exp(/3/ l (2r-7i)) (Bl) 

T = f> 

which can be performed without the explicit evaluation of the matrix elements H( p \(t,t). As we will see the trans- 
formed quantities will allow for a continuation to real n. The n — > limit then turns the finite dimensional eigenvalue 
equation into the integral equation (|49p. 

In section [TTT1 the entropic and the energetic parts of the Hessian, %s and He, were decomposed into a tensor 
product of identical 2x2 matrices. Therefore we first calculate the transformation (|B1[) for a tensor product of n 
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2x2 matrices g and use the general result for the two parts of the Hessian. 

n—p n 



E^HKg^h/ 5 ) exp (ph(2T-n)) 

r—p a—1 

n — p 

= E«+ih)'«+i )""'«- 



(B2) 



n 



a=l 
n — p 



g( 1+) I-) - h) 1+) Y \t P)„_ 2p exp ((3h(2r n)) 



» E«+ir"(H) B 



-|(-|g®g(l+) 



|g|+)(-|g|->-<+|g|-)(-|g|+)) E 

{Ta=±l} 
(T-p 



-cr-p 



n — 2p n — 2p 

E '-<-< 

ET a ^2r-n a=l 
n—c—p 



exp /3fe r a 



(T-p 



.6=1 



g+-g- 



g+-g- 



(+| g |r b )e^) Hgkc>e^< 

\r 6 =± / \t c =± 



c=l 

\ n — (7— p 



^^"fg + -^ + g__ e -^""" " 



where g CTT := (cr| g |t) , er, r G ± denote the matrix elements of g. In the second line we have used the definitions of 
{a; p\\ and |r; p) and then decomposed the symmetric part of |r; p) as in P2"| . using 2r — n = r a . 
The entropic contribution (|17p can be written as 



drdr' 



I exp {-ir'a 3 ) = J ■ 



2n c(r) / 27r c(r) 

v ' a=l J v ' a=l 

The entropic matrix g s is diagonal: g+ + = e~ %T , g? _ = e +4r , and using (|B2|) we find 



(B3) 



J2'H ( s p \a,T)e xp (l3h(2T-n)) = / _ ^ (^;p|| (g) g s (r') |r; p) exp (/fo(2r - n)) 

r— p "' r—p a—1 

= / ^^y(4 + ^)^ P (g S -e^)"-^ = / dr^*(r -(&,-»)) exp (^(2. - „)) 



E 



S<J,T 



c(2a — n) 



exp (/3h(2r - n)) . 



The expression in the last line leads to eq. P7|) . Due to the structure of the RS solution (Til)) it is convenient to 



determine the inverse of which amounts to 



The energetic contribution f|16p reads 

P n „ n 

H E = - dJGM{J)®[eP J &o+e-P J &i] =- dJG M (J)($ S E (J) 

J a=l a=l 

with diagonal elements g^ CT (J) = e' 3 ' 7 , and off-diagonal terms g^L ff = e _/3J 

n—p n—p n 

-^^ ) (a,r)exp(^(2r-n)) = j dJG M (J) E^HI^g^hrfexp^r-n)) 

J 1 6 

dJG Af (J) (d h u(h,J)) p w n (h,J) exp (/? (2cr - n) u(h, J)) 



= I dJG M (J)(e 2/3J ~e-^ J Y (jV+h) 



', a = ±1. Using (JB2J) we get 



-P(J+h)\ P ( p /3(-J+h) + e -P(-J+h)^ 



(B4) 



(B5) 



(B6) 
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with 

w n (h,J) = (2 cosh /3(J + ft)2 cosh /3(J - ft))' 

and 

1 tanh(^J) 



dhu(h, J) 



cosli (/3ft) 1 - tanlr (/3J) tanlr (^ft) 



Appendix C: The eigenvalue equation 

We start from the eigenvalue equation (|4"oj) for the matrices - H' p - ) . Each of the eigenvalues corresponds to one of the 
n + 1 — 2p representations Z?( p ) which arise in the decomposition (|37p . Splitting the eigenvalue equation into energetic 
and entropic part leads to 

n — p n — p 

A (P) X W (r) = J2 n { s p) (a, r) X (rt (r) + E «W (a, r) X ^ (a) (CI) 

<r—p a—p 

which is equivalent to 

d ~ h n~h) (p)f s _ , otfn „ ( p), , , /• dftpw 



A(P) / ^Sfer X (p) (r)exp(^ft(2r-n)) = X <'>(t) + / - E ^W) exp(/3ft^ - n)) X < P V) 

J (2cosh(pft)) n J (2 cosh(/?ft)) n 



Here we have used the explicit form of the inverse of tig given in (|B4I) . We now perform the transformations 
explained above which allows to perform the limit n — > 0. The transformation of the l.h.s. of the last equation reads 

r i "~ p r i 

A (p) /dftP(ft)- ^-^V exp(/3(ft + ft)(2r-n))x (p) (T) = / dft- ^-^V {h) <j> {p) {h + h) 

J (2cosh(/3ft)) ^ VV A >M w y (2cosh(/3ft)) K ' V ' 

whereas the r.h.s amounts to : 

n-jo . n-p 

Ex (p) Mexp(/3ft(2r-n))+ UhV(h) - E ^W) exp (/3(ft> ft)(2r - n)) X ( p V) 

r=p ^ (2cosh(/3ft)) 

/ n-p /n-p \ 
^W /o UafM* E E^W)exp(/?(ft + ft)(2r-n)) 
(2cosh(/3ft)) ^ p ^ y 

= (p) (ft) - JdhV(h) ( 2cQ 2(p~ h) y J dJG M (J) (d h u{h + h,J)Yw n (h + h,J)^ p) (u{h + h,J)), 
where we have used (|B6[) . The limit n — > may now be performed which yields our central result (|49l) . 
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